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Abstract 

In carrying out calculations of the effect of longitudinal space 
charge on longitudinal motion, the transverse beam size appears in 
a form factor which is usually written as (7 = 1 + 21n(5/a). In fact, 
this expression applies to particles with vanishing betatron amplitude 
in a beam with uniform transverse distribution. It is argued that an 
average over the transverse distribution should be used instead of the 
value on axis. It is shown that for the realistic ‘binomial’ family of 
distributions the 1 in the above expression for g should be replaced 
by a value near 0.5 if a is interpreted as twice the rms width of the 
beam. 


1 Introduction 

In an ideal machine with no space charge, the three degrees of freedom of 
motion can be investigated separately. With space charge, this is no longer 
possible: the longitudinal force on a particle depends upon its transverse 
position in the bunch and vice versa. It is still profitable to some extent to 
study longitudinal and transverse motions separately, but one must keep in 
mind that an approximation is being made. It has become standard practice 
to calculate longitudinal space charge effects using the on-axis force because 
it was felt that this corresponds to a ‘worst case’. However, there is no basis 
for this approach unless one is interested in the absolute longitudinal space 
charge limit (even in that case, the ‘on-axis’ approach is debatable). 
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Since synchrotron motion is slow compared with betatron motion, it is 
clear that in longitudinal simulations the longitudinal force should be aver¬ 
aged over transverse positions. (Conversely, in transverse simulations, the 
transverse force should be taken as that corresponding to an instantaneous 
longitudinal position.) This note is devoted to deriving a formula for the 
longitudinal force as averaged over transverse distributions for the binomial 
family of distributions. The beam is assumed to be round (radius=a) in a 
concentric round beam pipe (radius=6). It is found that the usual transverse 
geometric form factor [g), which is 1 -|- 21n(6/a) for the on-axis case of the 
uniform distribution, can be written as K + 2ln{b/a), where K depends upon 
the distribution and the dehnition of a. To be fair it should be pointed out 
that the error incurred by these assumptions (neither the beam nor the beam 
pipe is round nor of constant dimension) may very well be larger than the 
error incurred by using the on-axis formula. 


2 The Distribution 


The assumed distribution is 

p= + (1) 


for r < a and p = 0 otherwise. The normalization is // pdxdy = 2 ti f prdr = 
A, the line density. For p = 0 this is the Kapchinsky-Vladimirsky distribution 
(round case). For p —)■ cx) it is gaussian provided a —)■ cxd as well. This can 
be seen by hnding the standard deviation. 


-TT f pr^rdrdd 271 n , 

r2 = -—— = — / pr'^dr = -. 

/ prdrdO X Jo p + 2 


( 2 ) 


Since r'^ = x'^ + y'^ = x'^ + y'^ = 2x‘^ = 2al, this can be written as 


a = 


\/2(p -|- 2) ax- 


(3) 


In the Sacherer-Lapostolle convention, the beam size is characterized not by 
ax but by 5; = 2ax- This is handy because for the K-V case x = a. For the 
general case. 



( 4 ) 
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3 The On-Axis Potential 


From Gauss’ law, the electric field is radial with magnitude^ 

E{r) = i- /’■ pir'ydr' ^ ^ | 1 - (1 - f/aV' if r < a 

Cor JO 27reor [ 1 otherwi; 


( 5 ) 


With the beam pipe (radius= 6 ) at ground, the potential as a function of 
radius is 


G = 


E{r')dr'. 


( 6 ) 


Assuming all the beam is in the pipe {a < b), the potential in the beam is 


G = 


-A 


27ren 


'r[l-(l-rVa^) 




dr , (b 

— + In I - 

r 


a 


which for /i=integer can be expanded as 

-A 


V = 


47ren 




k=l 


k 


21 n r- 


(7) 


( 8 ) 


The restriction of jj, to an integer is of course not necessary, but gives simpler 
analytic expressions. 

The expression in square brackets can be identified as g. For the uniform 
distribution (/i= 0 ) on axis (r= 0 ), we recover 


g = 1 + 2 ln( 6 /a). (9) 

Generalizing to g = Kq + 21n(6/a), we find A'o = p See Table 1. 

For the gaussian, the sum diverges logarithmically, but for constant x, 
ln( 6 /a) diverges as well and the two infinities cancel.^ Defining Kq hy g = 
Ko + 2 \n{b/x), we get 

Kq = Kq — ln(/i /2 + 1 ). ( 10 ) 

Ko{+) is also given in Table 1. For large /i, Kq approaches 7 + ln(/r + 3/2), 
where 7 is Euler’s constant 0.57721566 ■ ■ ■. Hence, for the gaussian, Kq = 
7 + In 2 . 

^Strictly speaking this is only correct if A is constant. However it is a good approxi¬ 
mation if, as is usually the case with proton machines, the length scale of the longitudinal 
variation is large compared with the beam pipe size, i.e. bunch length 7> b. 

^Actually, this violates our assumption a < b. However, in realistic cases b > lax so 
there is so little beam outside r = b that the approximation is not significant. 
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Table 1: On-axis values for K 

II I ~ 


p 

Ko 

Ko 

0 

1.0000 

1.0000 

1 

1.5000 

1.0945 

2 

1.8333 

1.1402 

3 

2.0833 

1.1670 

5 

2.4500 

1.1972 

10 

3.0199 

1.2281 

00 

00 

1.2704 


4 The Average Potential 

One can find an average by first averaging the betatron motion for a given 
particle and then averaging over all the betatron amplitudes in the beam. 
However, for a stationary transverse distribution, one can find the result 
more easily by simply averaging the potential over all particles in the beam 
at a given instant in time without regard to betatron amplitudes. In this 
way the average potential is 


^ JVprdrde 
fprdrdO 


2 % /•“ 

T Jo 


V pr dr 


Integrating term by term we hnd 


(^) 


-A 

47reo 



k -\- p \ 


+ 2 In 


6 

a 


( 11 ) 


( 12 ) 


The summation can be identihed as {K). Again, it diverges as /i —)■ cxd, but 
remains hnite if the beam size is characterized by x instead of a. See Table 
2. For those interested in special functions, the solution for general p is 

(A) = 7 + 2^(/3) - ^(2/3 - 1 ) - ln(/3/2), (13) 


where 13 = p + 2 and ip is the psi (or digamma) function. An approximation 
that is within 0.2% for /i > 1 is (A) 7 — l/(4/i -|- 11). 


4 



Table 2: Average values for K 


/i 

(K) 

(K) 

0 

0.5000 

0.5000 

1 

0.9167 

0.5112 

2 

1.2167 

0.5235 

3 

1.4488 

0.5325 

5 

1.7968 

0.5440 

10 

2.3489 

0.5572 

oo 

oo 

0.5772 


5 Conclusion 

The value of K averaged over the beam is remarkably stable; varying by 
only 15% over all the binomial distributions from K-V to gaussian. It is 
therefore safe to conclude that for reasonably ‘good’ transverse distributions 
(i.e. smooth, peaked at the centre, no halo), the form factor g which appears 
in the longitudinal space charge force and in the space charge impedance is 
given to a good approximation by 

g = i + 2lu(l) (14) 

where b is the beam pipe radius and x is twice the rms beam size. This should 
be used in place oi g = 1 + 2 ln(6/ a) for the calculation of space charge effects 
when longitudinal motion is considered independently of transverse motion. 
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